Persistent currents in nanorings and quantum decoherence by Coulomb interaction 
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Employing instanton technique we evaluate equilibrium persistent current (PC) produced by a 
quantum particle moving in a periodic potential on a ring and interacting with a dissipative envi- 
ronment formed by diffusive electron gas. The model allows for detailed non-perturbative analysis 
of interaction effects and - depending on the system parameters - yields a rich structure of different 
regimes. We demonstrate that at low temperatures PC is exponentially suppressed at sufficiently 
large ring perimeters 27ri? > where the dephasing length is set by interactions and does 
not depend on temperature. This behavior represents a clear example of quantum decoherence by 
electron-electron interactions at T ^ 0. 



I. INTRODUCTION 

Electron decoherence is one of the key ingredients of 
the many-body ground state in the presence of disor- 
der and electron-electron interactions. The existing non- 
perturbative theory of this phenomenon at low tempera- 
tures in realistic disordered conductors^'^''^ is rather com- 
plicated, to a large extent because of the necessity to 
properly account for Fermi statistics for interacting elec- 
trons. At the same time the main physical reason for 
electron dephasing appears obvious already without un- 
necessary technical details: It is the electron interaction 
with the fluctuating quantum electromagnetic field pro- 
duced by other electrons moving in a disordered poten- 
tial. 

In order to be able to quantitatively describe and un- 
derstand the latter effect Guinea^ suggested a model 
which mimics all essential features of the "real" prob- 
lem of interacting electrons in a disordered conductor 
except for the Pauli exclusion principle. This model de- 
scribes a quantum particle moving on a ring with ra- 
dius R and interacting with quantum dissipative envi- 
ronment. For a system in thermodynamic equilibrium 
quantum decoherence manifests itself as effective sup- 
pression of off-diagonal density matrix elements beyond 
a certain length L^. Provided there exists nonzero elec- 
tron dephasing due to its interaction with quantum en- 
vironment at T — > 0, this dephasing length should 
stay finite down to zero temperature. Hence, all effects 
sensitive to quantum coherence, such as, e.g., persistent 
currents (PG) and Aharonov-Bohm (AB) oscillations in 
mesoscopic rings, should be suppressed by interactions as 
soon as the ring perimeter 27ri? exceeds L^. 

A great deal of information can be obtained by mod- 
elling the environment by a bath of Caldeira-Leggett 
(GL) oscillators. In this case it was demonstrated^ that 
PC is reduced by interactions in the ground state imply- 
ing suppression of quantum coherence exactly at T = 0. 
For the same GL environment Guinea^ found that AB os- 
cillations for a quantum particle on a ring are suppressed 
by the factor ~ exp{—{R/Lip)'^), where the length L^p is 
set by interactions and remains finite down to T = 0. A 
similar result was also obtained earlier from the real-time 



analysis^. Furthermore, the problem^i^ is exactly equiv- 
alent to that of Goulomb blockade in the so-called single 
electron box where exponential reduction of the effec- 
tive charging energy at large conductances^i^ is presently 
considered as a well established result. Thus, it is now 
widely accepted that PG for a quantum particle on a ring 
is exponentially reduced at large ring perimeters down to 
T = due to strong dephasing produced by interaction 
between the particle and the GL bath. 

It appears that presently no such consensus exists for 
another important model of the environment^ formed by 
a diffusive electron gas. Renormalization group argu- 
ments developed for this model^ suggest very weak power 
law ^ R^^ suppression of AB oscillations at T ^ 0, 
where the factor x 1 is set by interactions. On 
the contrary, the combination of semiclassics, instan- 
ton technique and quantum Monte Garlo (MG) analysis^ 
yields much stronger suppression of quantum coherence, 
namely exponential suppression exp{—R/Lp) (with 
temperature independent L^) at not too low T and power 
law suppression ^ R^^ with % w 1.8 at T — )■ and for 
2t:R > L^. 

More recently this problem was reconsidered by means 
of variational approach^*', perturbation theory^ ^ and MG 
simulations^^. In contrast to Ref. 9, either nc^*''^^ or 
very wcak^^ i?-dependent suppression of PG was found. 
Note, however, that the variational calculation^" reduces 
the PG problem to that of mass renormalization in the 
m ~ topological sector while, as we show here, PG 
is determined by other topological sectors (m ^ 0) that 
are distinct from the m = sector, unlike the varia- 
tional result in Ref. [lol . Perturbative in the interaction 
calculationsii can also miss the correct behavior of PG 
at not too small R as it was already demonstrated in 
Ref. for the problem under consideration and was also 
discussed elsewhere in a broader contextii^. More argu- 
ments along the same lines will be presented below in 
this paper. 

As far as numerical MG results are concerned, the 
authorsiS ascribed the difference between their conclu- 
sions and those of the previous work^ to insufficient Trot- 
ter number values employed in the MG analysis^. While 
this particular issue definitely requires further analysis, 
it is worth pointing out that the MG dataiS. cover only 
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FIG. 1: (Color online) The system under consideration: A 
particle on a ring in the presence of a periodic potential. The 
ring is pierced by the magnetic flux and the particle interacts 
with an effective environment formed by a dirty electron gas. 



the perturbative regime R < L^p where rather weak sup- 
pression of PC was found in Ref. |3 as well, cf. Fig. 3 in 
that paper. Hence, MC results^ for small R do not ap- 
pear conclusive for the most interesting non-perturbative 
regime R ^ where strong (though weaker than ex- 
ponential) interaction-induced suppression of PC was 
predicted- down to T ^ 0. 

Leaving detailed discussion of MC procedure and re- 
sults to the future, here we note that - despite signifi- 
cant theoretical activity in the field - currently no well- 
controlled non-perturbative approach is available which 
would allow to analytically study PC for the "particle -|- 
dirty electron gas" model- in the low temperature limit 
and at sufficiently large values of R. Given both funda- 
mental importance of the problem and remaining contro- 
versies in the literature, it is highly desirable to formulate 
such an approach which would unambiguously resolve the 
whole issue. A step in this direction is undertaken in the 
present work. 

The structure of our paper is as follows. In Sec. 2 
we define our model which is essentially identical to that 
pioneered by Guinea^ except we additionally introduce a 
periodic potential for the particle on a ring. In Sec. 3 
we formulate our main formalism and evaluate persistent 
current in the ring without interactions by means of the 
standard instanton technique. This approach is conve- 
niently generalized to the interacting regime in Sec. 4 
where we analyze suppression of both quantum coher- 
ence and PC by Coulomb interaction at sufficiently low 
temperatures down to T ^ and arbitrarily large ring 
perimeters. Discussion of our main observations and con- 
clusions is presented in Sec. 5. Some technical details of 
our calculation are specified in Appendices A and B. 



II. THE MODEL AND EFFECTIVE ACTION 

We will consider a quantum particle with mass M and 
electric charge e on a ring with radius R threaded by 
external magnetic flux see Fig. [T] As before^i^, it 
will be convenient to describe the particle position by 
a vector — {Rcos9, Rs'm9) and consider the angle 9 as 
a quantum variable. In contrast to Refs. HI^ Il0lllllll2l 
where the quantum particle on a ring was described only 
by its kinetic energy (i.e. no potential energy was in- 
cluded into consideration), in this paper we will assume 
that the particle moves in a periodic potential which - 
just for the sake of definiteness ~ is chosen in the form 
U{9) = Uo{l — cos{k9)). Here k is the total number of 
periods of the potential U{9) which the particle should 
pass before it makes one full circle on the ring. Accord- 
ingly, a non-interacting particle on a ring is described by 
the Hamiltonian 

Ho^^^^^^^^ + Uo{l-cos{n9)), (1) 

where $ = —i^od/d9 is the magnetic flux operator, 
Ec = l/(2Mi?^) and <i>o = 27rc/e is the flux quantum 
(here and below we set the Planck's constant equal to 
unity h = 1). 

Now let us include the interaction between the parti- 
cle on a ring and an effective dissipative environment. 
Specifically, we will assume that the ring is embedded in 
the environment formed by the so-called "dirty electron 
gas"—. The total Hamiltonian for our system reads 

H = Hq + Hci + Hint, (2) 

where H^i is the standard Hamiltonian for electrons in a 
disordered conductor and Hint describes interaction be- 
tween the particle and the electronic environment. Fluc- 
tuating electrons in this environment produce stochastic 
electromagnetic field V described by the equilibrium cor- 
relator 

{2nY k'^e{i\ujn\,k) 

where w„ = 2'KnT is the Matsubara frequency, e(u, k) is 
the dielectric susceptibility of the environment and X = 
R{t) - R{t'). Similarly to Refs. ii we will model the 
environment by 3d diffusive electron gas with 

1 ^ -iu + Dk"^ 
e{LO, k) Ana ' 

where a is the Drude conductivity of this gas, D — vpl/5 
is the electron diffusion coefficient and I is the electron 
elastic mean free path. Interaction between the particle 
on a ring and fluctuating electrons in the environment is 
described by the standard Coulomb term 

Hint = eV. (5) 
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In what follows we will assume that the whole system 
remains in thermodynamic equilibrium at a temperature 
T. Our first and standard step is to integrate out all envi- 
ronmental degrees of freedom effectively described by the 
collective variable V . In the limit of weakly disordered 
environment kpl ^ 1 to he analyzed below fluctuations 
of the field V can be considered Gaussian. In this case 
integration over this field is carried out exactly**'^. After 
that one arrives at the grand partition function of the 
system expressed as a single path integral over the angle 
variable 0{t): 

oc 
m— — oo 

= V / deal Veeyi^{2mm(j)^~ S[e]). (6) 



m— — oo 



Here we defined (3 = 1/T, cp^ = ^x/^o and the winding- 
number-projected partition functions 2^- The first term 
in the exponent in Eq. ^ takes care of the magnetic flux 
while the second term S[9] describes the effective action 
for our interacting particle on a ring. This action consists 
of two terms, 



S[0] = So[0] + Sint[0]- 



The term 



Jo 



1 fde 



AEc \dT 



+ C/o(l - cos(k6')) 



(7) 



(8) 



defines the action for a particle in the absence of the en- 
vironment. This action is identical to one for the Joseph- 
son junction^'^ where Ec plays the role of the charging 
energy. The term S'int describes the effect of interaction 
between the particle and the environment. For our model 
it has the form^i^ 

Sint[d\ =a dT dT . 2 , (9) 

Jo Jo sm [ttT[t-t')\ 

K{z) ^ 1 - ^ ^ (10) 



^Ar^ sin2(z/2) + l' 



where the constant a = 3/(8fc|.Z^) effectively controls 
the interaction strength in our model and r = R/l. Note 
that the integral in Eq. ^ is understood as a principal 
value. The formal divergence at r = r' is regularized 
by requiring K{Q) = which explains the origin of the 
first term in (fTO|) . For the sake of physical consistency of 
our model below we will set 1/kp I <^ 2ttR/k, where 
the first inequality just means that interaction should 
remain weak, a <^ I, while the second one implies that 
the distance between the neighboring potential minima 
should be much larger than I. Accordingly, the parameter 
r = R/l obeys the inequality 



r > k/2tt. 



(11) 



Note that it can also be convenient to rewrite the function 
K{z) in terms of the Fourier series 



K(z) — a„ sin^ 



nz 

T 



(12) 



where the Fourier coefficients are a„ ~ (2/7rr) ln(r/n) for 
1 < n r and a,i « otherwise. 



III. PERSISTENT CURRENT IN THE 
ABSENCE OF DISSIPATION 

Let us first evaluate persistent current in our system 
in the absence of interactions, i.e. we set the interaction 
constant equal to zero a = everywhere in this section. 
In what follows we will mainly be interested in sufficiently 
large values of the ring radius R. Hence, without loss of 
generality one can consider the limit Ec <C Uq. Below we 
will perform our calculation for a somewhat more strin- 
gent condition 



k'Ec < Uo 



(13) 



which simplifies our analysis but is by no means impor- 
tant for any of our key conclusions. 

In the tight binding limit (|13p and at sufficiently low 
temperatures the particle is located at the bottom of one 
of the potential wells (see Fig. [1]), i.e. in the vicinity of 
the points = 2'kp/k, where < p < k is an integer 
number. Accordingly, in Eq. ^ one should substitute 



2tt 



d0o 



E 



d0oS{0o - 2Trp/K) 



(14) 



The particle can move around the ring only by hopping 
between the neighboring minima = 2'itp/k, and = 
27r(p± 1)/k of the periodic potential U{0). Each of these 
tunneling events is described by the well known instanton 
(kink) trajectory 



0{t) — — arctan(exp(ajT)) 
and corresponds to the tunneling rate A/2, where 



(15) 



A = 



and uj = k\/2UqEc- In order to evaluate the grand par- 
tition function 



p=i 



-f3Ho 



(17) 



it is necessary to sum over all possible tunneling events 
of the particle between all potential minima to all or- 
ders in A. The minimum number of such hops should 
be equal to mn for any trajectory corresponding to the 
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winding number m. Taking into account that effective 



duration of each tunnehng event is 



and that the 



total imaginary time span equals to /3 = 1 /T we can dis- 
tinguish two different limits. In the limit tokw"^ <^ (3 the 
average distance between instantons is large as compared 
to their typical size, i.e. in this case we are dealing with 
dilute instanton gas. In the opposite limit mKw~^ 3> P 
instantons are very close to each other and essentially 
merge forming a single trajectory. Below we will also 
demonstrate that in the above conditions it suffices to 
set TO = 1. 



A. Dilute instanton gas 

We begin with the low temperature limit T <C w/k. To 
proceed let us evaluate a somewhat more general than in 
Eq. (|17p matrix element {pi\e~^^°\p2)- For this purpose 
we consider multi-instanton trajectories 



e(r) 



27rpi 



(18) 



where Vj = ±1 and Tj are respectively the topological 
charges and collective coordinates of instantons and 9{t) 
is defined in Eq. p5|) . The trajectory describes 
tunneling of the particle between the states \pi) and \p2) 
after to winds around the ring provided we fix 



n2 — ni — p2 — pi + rriK, 



(19) 



i.e. we consider configurations containing totally ni + n2 
instantons corresponding to ni hops clockwise and 712 
hops counterclockwise. Taking into account all possi- 
ble tunneling events restricted by the condition and 
summing over all winding numbers to, we obtain 



I — — 00 

(/5A/2) 



(20) 



m— — 00 
c>o /■ n A /r»\ni+n2 



E 



ni!n2! 



'^ni—n2,pi—p2- 



Making use of the integral representation for the Kro- 
necker symbol 



Zn 



(21) 



after performing a summation over m with the aid of 
Poisson's resummation formula 



E 



= E ^(^-^) 



(22) 



we obtain 



1/2 _^ 

e 2 



(23) 



This formula allows to easily recover the low-lying energy 
levels of our problem which contain all necessary infor- 
mation in order to evaluate PC. For instance, the ground 
state energy of the particle i?o('/'s) is obtained in a stan- 
dard way by taking the limit T = l//3— ^Oin Eq. 
which yields 



Eo{<l>x) - - - Acos ( — ^ 



(24) 



for -1/2 < (t)x < 1/2. Eq. ([24]) should be continued 
periodically outside this interval. This expression deter- 
mines the periodic dependence (with period equal to the 
flux quantum $0) of the ground state energy on the mag- 
netic flux ■ In the limit k 3> 1 this dependence reduces 
to a set of parabolas 

EQ{(t>x) - - ^^^^min„(02; - nf. (25) 

Turning back to the grand partition function (jl7p , from 
Eq. (HSl) we find 



(26) 



PC can now be easily obtained from the general formula 

eT d 



I = -- 



2tt 



■InZ, 



which yields diamagnetic current 
^_eA E;^iSin(g^) 



„/3Acos( 



J = l' 



/3Acos(- 



(27) 



(28) 



This expression fully determines PC in the ring at tem- 
peratures T <ti uj/ K and in the absence of interactions 
with dissipative environment. The dependence of the 
maximum PC value on temperature is depicted in Fig. 

m 

Eq. ((28|) is further simplified at temperatures above 
and below the interlevel distance ~ A/k^. In the limit 
T ^ A/k^ the leading contribution to the partition func- 
tion is defined by configurations with minimal number of 
instantons. Hence, in the sum over winding numbers 
in Eq. ^ it is sufficient to keep only the terms with 
TO — 0, ±1 terms. For the term with to = it is neces- 
sary to sum over all configurations, whereas for the case 
TO = ±1 only configurations with k instantons contribute. 
After some algebra we get 



/c— — 00 



Ate 



(29) 
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FIG. 2: (Color online) Maximum value of PC as a function 
of temperature for different values of k. Inset: The magnetic 
flux value (f)max at which PC reaches its maximum value as a 
function of T/A. 

where Io(a;) is the modified Bessel function of imaginary 
argument. For m = ±1 wc find 

= ^^e-'^-/^ (30) 

As a result we obtain 



with increasing T, cf. our Eq. ((32)) with Eq. (11) in 
Ref. dj. In the low temperature limit (T <C A/k^ here 
and T <C Ec for Uq — 0) the dependence I{(j)x) strongly 
deviates from sinusoidal (except for a special case k — 
In the limit k':^ 1 the sawtooth dependence ([34]) is iden- 
tical to that for the case C/q = where one has Ico ~ Ec- 
Finally, let us compare the dependence of the PC am- 
plitude /co on the ring radius R obtained in these two 
cases. Provided the parameter n is fixed and does not 
change with R, Eqs. and yield exponential de- 
cay of /co with increasing R since Ec oc 1/R^. This 
exponential decay is due to the fact that for larger R 
the potential profile changes in a way that the particle 
should tunnel at a longer distance 2'kR/ n between the 
two neighboring states \p) and |p ± 1). Alternatively, 
one can keep the distance between the adjacent poten- 
tial minima 6 — 2'iTp/ k and — 27r{p ± 1)/k unchanged 
while increasing R. This is achieved by varying the pa- 
rameter K with R as K (X R. Under this condition the 
tunneling rate A becomes independent of R and the 
magnitude of PC Ico (EH) decreases with increasing R 
as Ico oc 1/R^ exactly as in the case Uq = 0. Note, 
however, that due to the restriction p3p this regime can 
apply only at not too large values of R. 

B. Merged instantons 



/-/co(T) sin(2^0,). 



where 



IcoiT) 



i2T)-^K\lo{A/Ty 



(31) 



(32) 



At low temperatures T ^ A/k^ Eq. (pS)) reduces to a 
simple formula 

/=£^sinf^V -l/2<0,<l/2, (33) 

K \ K J 

which also trivially follows from Eq. ([24]) . This formula 
demonstrates that at T = the magnitude of PC in our 
system is proportional to A while its flux dependence 
deviates from the simple sinusoidal form for all k > 1 . In 
particular, in the limit k 3> 1 this dependence approaches 
a sawtooth one 



/ = /comin„( 



n), Ico = 27reA/K^ 



(34) 



Comparing the expressions for PC derived above with 
those for a free particle on a ring-i^iiSiiiii^ we observe 
that the main physical difference between these two mod- 
els is the presence of two distinct energy scales - to and A 
- in our case whereas only one energy scale Ec remains 
in the limit Uo ^ 0. Otherwise significant features of the 
effect are essentially the same in both models. Indeed, at 
high temperatures (T > A here and T > Ec for Uq = 0) 
the dependence of PC on the magnetic flux is sinu- 
soidal with the period $o and its amplitude decreases 



Now let us turn to the case of higher temperatures 
Lu/ K ^ T <^ Lu which can be realized in the limit k 3> 
1. In this case the leading contribution to the partition 
function originates from one multi-instanton trajectory. 
This trajectory of merged instantons 8^'"-'(r) can easily 
be evaluated due to the presence of the integral of motion 
which is the classical energy £m corresponding to the 
winding number m. This energy is fixed by the periodic 
boundary condition 



1 

T 



MR^ 



27r|m| 



de 



(35) 



and the trajectory 0(™)(r) is obtained from the equation 



T -To 



MR'^ 



de 



2 J ^£,„ + C/o(l - cos(k0)) 



(36) 



Making use of the standard quasiclassical technique one 
arrives at the following expression for the partition func- 
tion: 

where 



m— — oo 



d£l 



(37) 



S[£) = [ £ + Ua{l - cos{Ke))de (38) 

V Ec J 
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FIG. 3: (Color online) Different regimes for PC in the non- 
interacting case. The regions 1 (lowest T) and 11 (higher T) 
are described within the dilute instanton gas approximation 
(Eqs. (|33|l and (|3H) . p2|) respectively), whereas the high 
temperature region III corresponds to the case of merged in- 
stantons, Eqs. dSIJ, (gT]). The dashed lines T « A/k^ and 
T ~ uj/k indicate the crossover between these regimes. 



is the classical action. As before, the leading contribution 
to this partition functions comes from the lowest wind- 
ing numbers m = 0,±1. For m = one recovers the 
oscillator-like expression 



2 sinh 



(39) 



2T 



whereas for m = ±1 the instanton trajectory approaches 
the straight line in which case the integrals can be easily 
evaluated and yield 



Z±l = K 



2ttUoT _2xVot» 



(40) 



With the aid of these results one again arrives at the 
expression for PC in the form (j3ip with 



/co(r) = 2eT^— ^exp(^- j . (41) 

Different regimes considered in this section can also 
be summarized graphically by means of the diagram de- 
picted in Fig. O 



IV. EFFECT OF ELECTRON-ELECTRON 
INTERACTIONS 

A. Renormalization of the tunneling amplitude 

Now let us turn on interactions and analyze the effect 
of fluctuations in a dissipative environment. To this end 
we again employ the above instanton technique. Evalu- 
ating the path integral in ^ in the limit (fT51) and for 



T <^ uj/k we follow the same scheme and substitute the 
trajectory describing quantum tunneling of the par- 
ticle between different potential minima into the full ef- 
fective action ([7]). As before, let us fix the winding num- 
ber equal to m. Then we again arrive at configurations 
of totally k instantons restricted by the condition (fT5|) 
where we now also set pi = P2 = P- Evaluating the inter- 
action term 5int (©, (|10p on multi- instanton trajectories 
(fT8)l after some algebra (see Appendix A) we obtain 



S'int [©(''■)] = Airark/K 



-2a ^ VaVbgi'Pab) In 



a, 6—1, a<6 



sin[7rr(Tf, - To)] 



ttTuj 



(42) 



where 



g{(fi) = K{tp + 2tt/k) + K{ip - 2ti j k) - 2K{ip) (43) 



and 



27r ^ TT 

</5ab = > Vj (i/{, + Va) 

3=0. 



(44) 



We observe that the interaction term (1421) consists of two 
different contributions. One of them, Anark/n, describes 
interaction-induced suppression of quantum tunneling of 
the angle variable 6 between different potential minima. 
This term yields effective r-dependent renormalization of 
the tunneling amplitude 



Ae 



-47rQr/ K. 



(45) 



The remaining contribution in Eq. (j42p describes log- 
arithmic interaction between different instantons which 
occurs for k > 2 due to the presence of a dissipative envi- 
ronment. This logarithmic interaction is absent for k = 1 
in which case g{(p) = 0. 

At not too low temperatures T ^ Ar/n^ inter- 
instanton interactions just provide further renormal- 
ization of the tunneling amplitude Aj. — > Ar{l + 
2aK {2tt / k) ln{2TTT / uj)) . One can also write down the 
renormalization group (RG) equation 



dAr 
din a; 



2aK{2TT/K)Ar, 



(46) 



which yields both r— and T— dependent renormalized 
tunneling amplitude of the form 



A 



R 



2aK{2-!T/K) 



(47) 

Ar/K^ inter- 



At even lower temperatures T <^ 
actions again yield renormalization of the tun- 
neling amplitude which now becomes A^ — s- 

A,, (l + 2aK{2TT/K) In ^^J^) . The corresponding 

RG equation (|46p remains the same but relevant energy 
scale now becomes 



A„ = A. ^ 



2QJC(27r/^) 
l-2a/f (27r/K) 



(48) 
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Eq. (|48)) - together with Eq. p5|) - defines the renor- 
mahzed tunneling ampHtude in the Umit T — > and for 
K > 2. In the particular case k = 2 our results reduce to 
those established for the so-called spin-boson model with 
Ohmic dissipation^'''^^. 

We also note that Eq. (|48|) is formally applicable for 
a < 1/ {2K{2Tr/ k)), while for larger values of the interac- 
tion strength a and at T = the tunneling amplitude 
is renormalized to zero, A/j — 0. This is the conse- 
quence of the quantum dissipative phase transition which 
- similarly to other models— li^ii^ - occurs at the critical 
interaction strength ac = 1/{2K{2tt/k)) ~ 1/2 and im- 
plies localization of a quantum particle in one of potential 
wells at any a exceeding the critical value ac- Accord- 
ingly, no PC can flow at T = and a > 1 / {2K {27r / k)) . 
This formal conclusion, however, does not appear to be 
of substantial physical significance, since the applicabil- 
ity range of our model is restricted to small values of 
a < 1. 

Finally, we should point out that, employing a regu- 
lar perturbation theory in a, already in the first order 
one recovers additional terms which cannot be captured 
within the RG equation (I46p . The corresponding analysis 
is presented in Appendix B. 



B. Suppression of PC by interactions 

It turns out that the behavior of PC may be quite dif- 
ferent depending both on temperature and on the param- 
eter K. Therefore, it is appropriate to distinguish several 
different cases. 



1. One potential minimum k — 1 

As we already discussed, in the case k = 1 logarithmic 
interaction between instantons is absent, and the only 
effect of interactions is r— dependent renormalization of 
the tunneling amplitude (I45p . Accordingly, for PC in this 
case we obtain 



/ eAe 



' sin(27r(/)^). 



(49) 



This result is valid at all temperatures in the range 
T <C w. It demonstrates that for k = 1 Coulomb in- 
teraction yields exponential suppression of PC down to 
T = provided the ring perimeter 2t:R exceeds an effec- 
tive dephasing length 



which is set by the effective interaction strength a and 
does not depend on temperature. Note that exactly the 
same length scale was found in the absence of the periodic 
potential [/o in Ref. '9'. 

2. Not too low temperatures and k > 2 



In the case k > 2 instantons interact logarithmically 
and - in addition to (|45p - at not too low temperatures 
the tunneling amplitude is renormalized according to Eq. 
p7)l . Combining these two equations and substituting 
the renormalized tunneling amplitude An into Eq. 
instead of A we obtain 



where 



/c(r) sin(27r0,). 



(51) 



(50) 



^ ^ I /co(r)e-4-"-, T^u/k 

(52) 

and /co {T) in the corresponding limits is defined respec- 
tively in Eq. p2|) and ((4T|) . In this case we again re- 
cover the same exponential suppression of PC at ring 
perimeters 2'kR ^ L^p with temperature independent de- 
phasing length defined in Eq. (|50p . At the same time, 
the pre-exponent in the expression for Ic depends on 
temperature as a power law Ic{T) oc T"^ with fi ~ 
k(1 — 2aK {2tt / k)) ~ 1. For k > 2 and small values of 
the interaction strength a <C 1 we have n > 0, i.e. Ic{T) 
grows with decreasing temperature. This growth, though 
somewhat weaker than in the non-interacting case (|32p 
(since fi < k — 1), implies that Eqs. (|5T|) . ([52]) can be 
trusted only at T > A^j/k^ whereas at even lower tem- 
peratures one expects a crossover to a different regime to 
be discussed below. We also note that qualitatively sim- 
ilar behavior of PC at not too low temperatures follows 
from the numerical analysis^ of the model with Uq = 0. 



3. Zero temperature limit and k > 2 

In order to evaluate PC in the limit T ^ A^/ we will 
make use of Eq. (jB14[l for the free energy. Combining 
this expression with Eq. (|27p and observing that the 
difference A^ — A^ ~ a <C 1 in the last term in Eq. 
(|B14p can be safely neglected within the accuracy of our 
calculation (since it only produces extra terms ~ a^), we 
obtain 



2nQ 



sm 



n— — r 



a„ sm [^—j 



cos 



/ 2'K(l>x 



In 



s,n (-J sm 



(53) 
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where A/j is defined in Eqs. (|48|) . (|45p. This result allows 
to make the following observations. 

Firstly, taking into account the dependence of A/j on 
r we conclude that exactly at T = and at sufhciently 
large ring perimeters PC is exponentially suppressed as 

i.e. in this case for any fixed value k we can define an 
effective zero temperature dephasing length 

= - 2aX(27r/K)). (55) 

This result demonstrates that in the limit of weak inter- 
actions a ^ 1 the effective length (|55|) turns out to be 
approximately k times longer than L^p^ i.e. L^p « l/an. 
Should, however, we assume that k (x r, no finite dephas- 
ing length could be defined from Eq. (|54p . although also 
in this case even for small a PC can suffer exponentially 
strong suppression by electron-electron interactions due 
to the condition (fTT|) . 

Secondly, the result (|53p demonstrates that at T ^ 
the effect of electron-electron interactions on PC does not 
just reduce to renormalization of the tunneling amplitude 
A — * Aij. We observe that Eq. also contains addi- 
tional terms evaluated here within the first order pertur- 
bation theory in a. This r-dependent first order contribu- 
tion turns out to be singular at values of (jix close to ±1/2 



and at all other half-integer numbers which indicates in- 
sufficiency of the first order perturbation theory, at least 
for such values of (px- Higher order terms of the per- 
turbation theory in the interaction may contain similar 
(or even stronger) singularities and, on top of that, may 
grow with increasing r. Hence, at T — > no perturbation 
theory in a can in general be trusted, in particular at suf- 
ficiently large r. A more detailed analysis of this issue 
is beyond the scope of the present paper. Here we only 
conjecture that such analysis might yield an additional 
dependence of PC on the ring radius r not accounted for 
in the expression for A/j. It is quite likely that such kind 
of r-dependence of PC at T = was also observed within 
a numerical treatment^ in the case Uq — Q. Additional 
support for this conjecture is provided by the exact solu- 
tion presented below. 



4.. Toulouse limit 

As we already pointed out, for k — 2 our problem is 
exactly mapped onto the well known spin-boson model 
with Ohmic dissipatio n^^i^^ . In this case interaction be- 
tween instantons is given by Eq. (|A9|) and the grand 
partition function reads (see also Appendices A and B) 
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Z[z;p] = (Arcosz)^" / dn / dra... / dT2„exp 4a ^ {-ly+Hn 

n=0 i i . \ :)<k=l 



(56) 



r 



In order to non-perturbatively evaluate PC at all values 
of r and at all temperatures including T = we can profit 
from the exact solution known for the particular value of 
the interaction strength a = 1/4, the so-called Toulouse 
limit. 



Introducing the parameter 



1 



cos z 



one can 



conveniently write down the exact expression for the par- 
tition function (j56|) with a = 1/4 in the form 



Z[z; l3] — exp 



dx- 



1-e 2^e<^) 



: sinh ; 



7rT/i^ 



This result allows to immediately establish PC which 
reads 



I = ^—^ sin(27r(/)^) dx- 

Zuj / smh X 



^ cos (ir0^) 



(58) 



At not too low temperatures lj ^ T ^ A^ we obtain 

2ijj 



eA^ 

I = e-^'' In sin(27r0^). 

2u) ttT 



In low temperature limit T one finds 



I = -^sm{2TT(t)x) I dy- 



(59) 



(60) 



(57) 

Evaluating the integral in (j60p we arrive at the final result 



eA2 
~2c7 



,2.-7 



ttA^ cos^ (Tn/)^, ) 



7rr + In ( — — — — ; — r 1 sin(27r(/)a;). 

(61) 

We observe that both at non-zero temperatures and ex- 
actly at r = the above exact expressions for PC demon- 
strate its exponential suppression at ring perimeters ex- 
ceeding = (these two length scales coincide for 
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K = 2 and a = 1/4, see Eq. ([55]) ). In addition, the re- 
sult (pT|) demonstrates that at T = (a) the dependence 
of PC on r deviates from purely exponential (which is in 
agreement with our above conjecture) and (b) the depen- 
dence of PC on deviates from purely sinusoidal and 
contains the logarithmic singularity at half-integer values 
of cf. also Eq. ([55)1 . 

V. DISCUSSION 

In this paper we proposed a model which allows 
for detailed non-perturbative treatment of the effect of 
electron-electron interactions on PC in normal nanorings 
at low temperatures. Our investigation employs a well- 
controlled instanton technique and yields a rich struc- 
ture of different regimes. The main features observed 
within our analysis can be summarized as follows: (i) 
Coulomb interaction yields i?-dependent renormalization 
of the tunneling amplitude (I45p which, in turn, results in 
exponential suppression of PC at large enough _R, (ii) log- 
arithmic interaction between instantons yields additional 
renormalization of the tunneling amplitude described by 
the RG equation (|^^ and {Hi) electron-electron interac- 
tions generate yet additional contributions not captured 
by Eqs. (gS]) and (ge]), see Eq. ((531) and Appendix B. 
These contributions may become particularly important 
at T ^ indicating the failure of the naive perturbation 
theory in the interaction at sufhciently large R. 

Although the effect (Hi) still requires additional non- 
perturbative analysis, already (i) and (ii) result in expo- 
nential suppression of PC at any temperature including 
T —f for any given k and at ring perimeters exceeding 
the dephasing length set by interactions and defined in 
Eqs. (|50p and ([55]) . Note that the length scale identi- 
cal to ((50)) also follows from the earlier non-perturbative 
analysis^, developed for the case Uq — 0. Thus, similarly 
to Ref. [§ the decoherence effect in our model is controlled 
by the parameter 

r 

n=l 

rather than by a or a In r as it was sometimes suggested 
in the literature in the case Uq ~ 0. 

It is worthwhile to stress that exponential dependence 
of PC on R of the form / oc exp{—AR) by itself does not 
yet necessarily imply decoherence. For instance, even 
in the absence of interactions PC / oc A ([33]) can de- 
crease exponentially with increasing R provided the pa- 
rameter K is fixed to be independent of the ring radius. 
Obviously, quantum coherence is not destroyed in this 
case. Exponential reduction of PC with increasing R at 
T ^ can also occur in superconducting nanorings due 
to proliferation of quantum phase slips^^ii^. Also in that 
case the dependence / oc exp(— yli?) can be interpreted 
just as a non-trivial coordinate-dependent renormaliza- 
tion effectJi. 



An important qualitative difference between nanorings 
with dissipation considered in Sec. 4 and the two last ex- 
amples is that in our problem dissipation explicitly vio- 
lates time-reversal symmetry (thus causing genuine deco- 
herence of a quantum particle), while no such symmetry 
is violated in superconducting nanoring a^^'^^ or in the 
absence of dissipation (Sec. 3). Hence, quantum coher- 
ence remains fully preserved in the last two cases despite 
exponential suppression of PC at large R. 

One can discriminate between decoherence and pure 
renormalization in a number of ways. For instance, one 
can drive the system out of equilibrium and investigate 
its relaxation by means of a real-time analysis. This ap- 
proach was employed in Refs. [gllTsI where a finite de- 
phasing time Tip = Lip/v was found at T = with 
defined in Eq. ([50)) and v being the particle velocity. 
This observation allows to unambiguously identify quan- 
tum decoherence. 

Another wayS- amounts to analyzing fluctuations of PC 
in the ground state of an interacting system. For in- 
stance, one can study the correlator ((/ — (/))^), where 
/ is the current operator which expectation value (/) de- 
fines PC in the ground state. Within the model^ it was 
demonstrated that, while the average PC in the ground 
state decreases with increasing interaction strength, its 
fluctuations increase, thus implying genuine decoherence 
rather than pure renormalization. A similar situation 
occurs within the model studied here. 

Following^ here we intentionally disregarded Fermi 
statistics by suppressing electron exchanges between the 
ring and the environment. The question arises if inclusion 
of the Pauli principle into the model could alter our main 
conclusion about non-vanishing electron decoherence at 
zero temperature. Golubev and one of the present au- 
thors addressed this issue by developing two entirely dif- 
ferent non-perturbative in the interaction techniquesii^. 
Both these methods yield the same conclusion: Although 
Fermi statistics is crucially important for other proper- 
ties of a dirty interacting electron gas, it practically does 
not affect interaction-induced quantum decoherence at 
T — !■ 0. Quantitative agreement was demonstrated^- be- 
tween the corresponding results derived for the models 
within and without^ the Pauli principle. 

Despite all these developments (including recently ob- 
tained exact solution of the problem^) it is sometimes 
argued in the literature that the Pauli principle can pre- 
clude from electron decoherence at T ^ 0. E.g. this con- 
clusion was reached on the basis of the first order pertur- 
bation theory in the interaction^^. Insufficiency of such 
kind of perturbation theory for the problem in question 
was repeatedly demonstrated elsewhereii^'^'^^ and was 
also observed here within the model considered. More 
recently, von Delft and co-authors^^ re-iterated the same 
incorrect conclusioni^. Unfortunately the analysis^! is 
fundamentally flawed since it violates causality effectively 
implying that dynamics of interacting electrons would be 
affected by photons coming both from the past and from 
the future. Hence, the results^ cannot be considered se- 
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APPENDIX A: INSTANTON GAS IN THE PRESENCE OF INTERACTIONS 



Let us analyze the effect of electron-electron interactions on the dilute instanton gas employed in our work. Com- 
bining Eqs. ((HI) and (|12p we can rewrite the dissipative part of the action describing such interactions within our 
model in the following form: 

(3/2 (3/2 

S:nt[e]^--Y.an / dT / dr' ■ 2, J (^1) 

^ J J sm (ttJ (t — r)) 

-13/2 -[3/2 

Substituting multi- instanton trajectories 0(t) (jlSp into Eq. (|Aip and employing the approximation 

1 - e"®(^) « e~ £ (1 _ ^^nui9(r-r,)•^ ^ ^ Cj/„j(t) (A2) 

'j 3 

appropriate for well-separated instantons, with the aid of the identity 

-^T^t ^ - 1^ sin(^r(r - r'))| (A3) 

sm (ttJ (r — t')) otot 

we obtain 

13/2 13/2 

5int[e] = -|^a„^C,,C,, j dr j dT'/;,,(r)4.^(r')ln|sin(^r(r-r'))|. (A4) 

i^'i^ -0/2 -S3/2 

Let us first consider the terms in the above sum which correspond to j\ = 22 ■ They are 

/9/2 /3/2 

-fE«"/ '^^ I dr'/;,(r)4-(r')ln|sin(^r(r~r'))|. (A5) 

-13/2 -13/2 

Since an effective instanton width 1/lu is much smaller than the inverse temperature /3 it is possible to split the 
contribution (jA5[) into regular and singular (in the limit T — ^ 0) parts, respectively 



00 oc 



5.cg = - f E / / '^^'aS^ cos(n(0~(r) - dir'))) In \u;{t - r')\ ^ ^ (A6) 



-00 —00 



and 



Q . — _ 



2a > a„ sin^ ( — In = -2aK — In — . (A7) 

^-^ \ K / UJ \ ^ J ^ 



Consider now the remaining terms with ji 7^ j2- Provided the distance between instantons remains much larger than 
their width ~ l/w we find 

- «E«" E ^ [q,Q,(1 - e^^^O(l - e-^''^0] ln(sin(^T(r,, - t,J)). (A8) 

" 'h<h 
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Performing summation over n in Eq. (jASp and combining the result with Eqs. (|A6|) and (IA7|) we arrive at Eq. (|42)) . 

Let us note that in the particular case k — I the logarithmic interaction between instantons vanish, while in the 
case of K = 2 and for large enough r we get ^intl©] = 27rar(ni + ^2) + S{a) where the logarithmic inter-instanton 
interaction S{a) takes the form 

Sia) = -4a ^ (-1)-- m ^in(.T(r,._-.,J N _ ^^^^ 

Hence, the partition function for our model with k = 2 is formally equivalent to that for the well known spin-boson 
model with Ohmic dissipationi^i^. 



APPENDIX B: PERTURBATION THEORY 



Let us rewrite the partition function of our problem as 



g27rim</)^-Si„t[6(T)] J. 



Ti r„_i 

With the aid of the Poisson's resummation formulae Eq. (|Bip can be transformed to 

2ir „ 



/d7 
n m=— 00 k=l 



where 



^[-;/5]-E E ■• E 



Ti T„ _ 1 



and 



W^(Ta - n) = In 



/3 ;3 /3 fc fc 

dri / dT2... / dr„e 



sin[7rr(Tfc - Tg)] 



n— — oo 



(Bl) 



(B2) 



(B3) 



(B4) 



Expanding the exponent in Eq. (IB3p one recovers the perturbation series for Z[z] [3]. Extending the definition of the 
Fourier coefficients a„ to negative n in such a way that a_„ = a„ and ao = we may write 



g{ipab) = 2 ^ a„ sin^ cos{nipab) = ^ a„ sin^ (J^'j ^-^n^at^ 



(B5) 



With the aid of this equation it is easy to define all orders of the perturbative expansion Z — Z^"^ + Z^^^ + where 

ZW[z;/3] =e-2^^^'''"'(^/2) (B6) 

and 

P 



Z^^'^[z-(3] = -2aAl a„sin2 (2 + ^) sin^ {^'^ J dn J dT2Zo[z:Ti]Zo 



27171 

z H ;t2 - Ti 

K 



Za[z;[3-T2]W{T2-Tl). 

(B7) 

The last term can also be represented graphically by the diagram in Fig. [H Performing the Fourier transformation 
of the partition function 



d\ 



Z[z,(i]^ / — e-^/3^[z,A] 



27r 



(B8) 
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FIG. 4: (Color online) First order diagram of the perturbation theory in a. 



and introducing the self-energy E as a sum of all irreducible diagrams, we obtain 

Z[z, A] ^ 



-i\ + 2A^ fAn\z/2) - S[z; A] 



To the first order in the interaction one finds 











sm^ 






K / 





k — — OG 



iX - 2TrikT + 2Ar sin^(§ + ^) 



(B9) 



(BIO) 



where Wq — ln[27rT/a;] and Wk — l/2\k\ for k ^ 0. Performing summation over fc, we get 
El [z; A] = 2aAl ^ o„ sin^ (^z + — ^ sin^ ^ 

n——r 

In If + 7 + ^V^ (l + ^+2'A--"'(t+^) 



2ttT 



i^(l- 



27rT 



(Bll) 

-?;A + 2A^sin^(f + ^ ^ 

where ip is the digamma function and 7 ~ 0.577 is the Euler constant. Substituting this expression into Eq. (IB9p . we 
observe that Z[z, A] has a pole at Ap = — 2iAr sin^(z/2) + 6Xp, where 



SXp — iaAr a„ sin 





■Kn\ 












X 











In 



2ttT 



+ 7+ -V' 1 + 



■ sm z 



7rn\ 



sm 



TTJl 



+ 2^(1 



iAr 



■ sm 



/ 7m\ /"!m\ 



(B12) 



Combining the above results with Eq. (|B8p one arrives at the final perturbative in the interaction expression for the 
partition function. Of particular interest is the low temperature limit T <C A^/k^. In this case one has '0(1 + z) ~ 
ln(z) + 0{l/z) and the free energy reduces to 

27r0T. + TT 



TlnZ = 2ArSin I \ — aAr / a„ sin ^ — J sir 



, 2Ae'^ , 
In h In 



sm I — J sm 



2nQ 



(B13) 



It is straightforward to observe that the first logarithmic term in the square brackets simply yields the renormalization 

A,- Ar ^1 + 2aK{2TT/ k) In ^^J^^ and does not alter the dependence of PC on the flux (p^. Hence, this perturbative 

contribution can be absorbed in the first term in (jB13p simply by substituting A/? instead of A^, where Ajj is defined 
in Eq. (148]) . 

The second logarithmic term in the square brackets in (jBlSP contains an additional flux dependence which turns 
singular at (j)^ close to the half-integer numbers. Clearly, this perturbative contribution cannot be just reduced to the 
renormalization (I48p . As a result, we obtain 



- TlnZ = 2Ai^sin2 f ^ j - aAr ^ a„ 

\ / n——r 

Employing this expression we arrive at Eq. (|53p . 



, sm 



V K 



sm 



27rfi 



In 



sm 



sm 



27rfi 



(B14) 
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